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RFFT gFNT TMAGF. RRCQNSTRUCTION ALGORITHM FOR JHE 
CIRCLE AND LINE CONE BEAM COMPUTED TOMOGRAPHY 



This invention is a Continuation-In-Part of United States Patent Application Serial 
5 No. 1 0/728, 1 36, filed December 4, 2003 which claims the benefit of priority to U.S. 
Provisional Application, Serial No. 60/430,802 filed, December 4, 2002, and is a 
Continuation-In-Part of United States Patent Application Serial No. 10/389,534 filed 
March 14, 2003 which is a Continuation-In-Part of Serial No. 10/389,090 filed March 14, 
2003, which is a Continuation-In-Part of Serial No. 10/143,160 filed May 10, 2002 now 
10 U.S. Patent 6,574,299, which claims the benefit of priority to U.S. Provisional 
Application 60/312,827 filed August 16, 2001. 

FIELD OF INVENTION 
This invention relates to computer tomography, and in particular to processes, 
1 5 methods and systems for reconstructing three dimensional images from the data obtained 
by a circle and line scan of an object, such as when the object first moves through a C- 
arm, and then the C-arm rotates around either a non moving or a moving object. 

BACKGROUND AND PRIOR ART 
20 Over the last thirty years, computer tomography (CT) has gone from image 

reconstruction based on scanning in a slice-by-slice process to spiral scanning to also 
include non-spiral scaiming techniques such as those performed with C-arm devices, with 
all techniques and devices experiencing problems with image reconstruction. 
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From the 1970s to 1980s the slice-by-slice scanning was used. In this mode the 
incremental motions of the patient on the table through the gantry and the gantry rotations 
were performed one after another. Since the patient was stationary during the gantry 
rotations, the trajectory of the x-ray source around the patient was circular. Pre-selected 
5 slices through the patient were reconstructed using the data obtained by such circular 
scans. 

From the mid 1980s to present day, spiral type scanning has become the preferred 
process for data collection in CT. Under spiral scanning a table with the patient 
continuously moves at a constant speed through the gantry that is continuously rotating 

10 about the table. At first, spiral scanning has used one-dimensional detectors, which 
receive data in one dimension (a single row of detectors). Later, two-dimensional 
detectors, where multiple rows (two or more rows) of detectors sit next to one another, 
have been introduced. In CT there have been significant problems for image 
reconstruction especially for two-dimensional detectors. Data provided by the two- 

15 dimensional detectors will be referred to as cone-beam (CB) data or CB projections. 

In addition to spiral scans there are non-spiral scans, in which the trajectory of the 
x-ray source is different from spiral. In medical imaging, non-spiral scans are frequently 
performed using a C-arm device, which is usually smaller and more portable than spiral 
type scanning systems. For example, C-arm scanning devices have been useful for being 

20 moved in and out of operating rooms, and the like. 

Fig. 1 shows a typical prior art arrangement of a patient on a table that moves 
through a C-arm device, that is capable of rotating around the patient, having an x-ray 
tube source and a detector array, where cone beam projection data sets are received by the 
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x-ray detector, and an image reconstruction process takes place in a computer with a 
display for the reconstructed image. 

There are known problems with using C-arm devices to reconstruct data. See in 
particular for example, pages 755-759 of Kudo, Hiroyuki et al., Fast and Stable Cone- 
5 Beam Filtered Backprojection Method for Non-planar Orbits, lOP Publishing LTD, 1 998, 
pages 747-760. The Kudo paper describes image reconstruction using C-arm devices for 
various shift-variant filtered back projection (FBP) structures, which are less efficient 
than convolution-based FBP algorithms. 

For three-dimensional (also known as volumetric) image reconstruction from the 
10 data provided by spiral and non-spiral scans with two-dimensional detectors, there are 
^ two known groups of algorithms: Exact algorithms and Approximate algorithms, that 

each have known problems. Under ideal circumstances, exact algorithms can provide a 
replication of an exact image. Thus, one should expect that exact algorithms would 
produce images of good quality even under non-ideal (that is, realistic) circumstances. 
1 5 However, exact algorithms can be known to take many hours to provide an image 

reconstruction, and can take up great amounts of computer power when being used. 
These algorithms can require keeping considerable amounts of cone beam projections in 
\ memory. 

Approximate algorithms possess a filtered back projection (FBP) structure, so 
20 they can produce an image very efficiently and using less computing power than Exact 
algorithms. However, even under the ideal circumstances these algorithms produce an 
approximate image that may be similar to but still different from the exact image. In 
particular. Approximate algorithms can create artifacts, which are false features in an 
image. Under certain circumstances and conditions these artifacts could be quite severe. 
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To date, there are no known algorithms that can combine the beneficial attributes 
of Exact and Approximate algorithms into a single algorithm that is capable of replicating 
an exact image under the ideal circumstances, uses small amounts of computer power, 
and reconstructs the exact images in an efficient manner (i.e., using the FBP structure) in 
5 the cases of complete circle and line and incomplete circle and line scanning. 

If the C-arm rotates 360 degrees around the patient, this produces a complete 
circle. If the C-arm rotates less than 360 degrees around the patient, this produces an 
incomplete circle. In what follows, the word circle covers both complete and incomplete 
cases. Here and everywhere below by the phrase that the algorithm of the invention 

10 reconstructs an exact image we will mean that the algorithm is capable of reconstructing 
an exact image. Since in real life any data contains noise and other imperfections, no 
algorithm is capable of reconstructing an exact image. 

Image reconstruction has been proposed in many U.S. Patents. See for example, 
U.S. Patents: 5,663,995 and 5,706,325 and 5,784,481 and 6,014,419 to Hu; 5,881,123 and 

15 5,926,521 and 6,130,930 and 6,233,303 to Tarn; 5,960,055 to Samaresekera et al.; 
5,995,580 to Schaller; 6,009,142 to Sauer; 6,072,851 to Sivers; 6,173,032 to Besson; 
6,198,789 to Dafni; 6,215,841 and 6,266,388 to Hsieh. Other U.S. patents have also been 
proposed for image reconstruction as well. See U.S. Patents: 6,504,892 to Ning; 
6,148,056 to Lin; 5,784,481 to Hu; 5,706,325 to Hu; and 5,170,439 to Zeng et al. 

20 However, none of the patents overcome all of the deficiencies to image 

reconstruction referenced above. The inventor is not aware of any known processes, 
methods and systems that combines the beneficial attributes of Exact and Approximate 
algorithms into a single algorithm that is capable of replicating an exact image under the 
ideal circumstances, uses small amounts of computer power, and reconstructs the exact 
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images in an efficient manner (i.e., using the FBP structure) in the cases of complete 
circle and line and incomplete circle and line scanning. 

SUMMARY OF THE INVENTION 
5 A primary objective of the invention is to provide improved processes, methods 

and systems for reconstructing images of objects that have been scanned with two- 
dimensional detectors. 

A secondary objective of the invention is to provide improved processes, methods 
and systems for reconstructing images of objects scarmed with a circle and line x-ray 
10 source trajectory that is able to reconstruct an exact image and not an approximate image. 

A third objective of the invention is to provide improved processes, methods and 
systems for reconstructing images of objects scanned with a circle and line x-ray source 
trajectory that creates an exact image in an efficient manner using a filtered back 
projection (FBP) structure. 
15 A fourth objective of the invention is to provide improved processes, methods and 

systems for reconstructing images of objects scanned with a circle and line x-ray source 
trajectory that creates an exact image with minimal computer power. 

A fifth objective of the invention is to provide improved processes, methods and 
systems for reconstructing images of objects scanned with a circle and line x-ray source 
20 trajectory that creates an exact image utilizing a convolution-based FBP structure. 

A sixth objective of the invention is to provide improved processes, methods and 
systems for reconstructing images of objects scanned with a circle and line x-ray source 
trajectory that is CB projection driven allowing for the algorithm to work simultaneously 
with the CB data acquisition. 
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A seventh objective of the invention is to provide improved processes, methods 
and systems for reconstructing images of objects scanned with a circle and line x-ray 
source trajectory that does not require storing numerous CB projections in computer 
memory. 

5 An eighth objective of the invention is to provide improved processes, methods 

and systems for reconstructing images of objects scanned with a circle and line x-ray 
source trajectory that allows for almost real time imaging to occur where images are 
displayed as soon as a slice measurement is completed. 

A preferred embodiment of the invention uses a seven overall step process for 

10 reconstructing the image of an object under a circle and line scan. In a first step a current 
CB projection is measured. Next, a family of lines is identified on a detector according to 
a novel algorithm. Next, a computation of derivatives between neighboring projections 
occurs and is followed by a convolution of the derivatives with a filter along lines firom 
the selected family of lines. Next, using the filtered data, the image is updated by 

15 performing back projection. Finally, the preceding steps are repeated for each CB 

projection until an entire object has been scanned. This embodiment works with keeping 
several (approximately 2-4) CB projections in memory at a time and uses one family of 
lines. 

The invention is not limited to an object that undergoes a scan consisting of a 
20 single circle and a single line segment. The invention can be applied to trajectories 

consisting of several circles and line segments by applying it to various circle and line 
pairs, and then combining the results. 

The circle and line scanning can include a partial planar curve scan before or after 
a line scan. The planar curved scan can be less than a full circle and even greater than a 
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full circle. Additional and subsequent circle and line scans can be done consecutively 
after a first circle and line scan. 

Further objects and advantages of this invention will be apparent from the 
following detailed description of the presently preferred embodiments, which are 
5 illustrated schematically in the accompanying drawings. 



BRIEF DESCRIPTION OF THE FIGURES 
Fig. 1 shows a typical prior art view arrangement of a patient on a table that moves 
through a C-arm device, that is capable of rotating around the patient, having an x-ray 
10 tube source and a detector array, where cone beam projection data sets are received by the 
x-ray detector, and an image reconstruction process takes place in a computer with a 
display for the reconstructed image. 

Fig. 2 shows an overview of the basic process steps of the invention. 
Fig. 3 shows mathematical notations of the circle and line scan. 
15 Fig. 4 illustrates a stereographic projection from the current source position on to the 
detector plane used in the algorithm for the invention. 

Fig. 5 illustrates finding of a filtering line for a reconstruction point x when the x-ray 
source is on the line. 

Fig. 6 illustrates finding of a filtering line for a reconstruction point x when the x-ray 
20 source is on the circle. 

Fig. 7 illustrates a family of lines used in the algorithm of the invention corresponding to 
the case when the x-ray source is on the line. 

Fig. 8 illustrates a family of lines used in the algorithm of the invention corresponding to 
the case when the x-ray source is on the circle. 



. p ^ ..J... y [j J jg Q 114. ;g 5; 3 s p^^y Y)]^^ . 

UCF-456CIP 

Fig. 9 is a three substep flow chart for identifying the set of lines, which corresponds to 
step 20 of Fig. 2. 

Fig. 10 is a seven substep flow chart for preparation for fihering, which corresponds to 
step 40 of Fig. 2. 

5 Fig. 1 1 is a seven substep flow chart for fiUering, which corresponds to step 50 of Fig. 2. 
Fig. 12 is an eight substep flow chart for backprojection, which corresponds to step 60 of 
Fig. 2. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
10 Before explaining the disclosed embodiments of the present invention in detail it 

is to be understood that the invention is not limited in its appHcation to the details of the 
particular arrangements shown since the invention is capable of other embodiments. 
Also, the terminology used herein is for the purpose of description and not of limitation. 

This invention is a Continuation-In-Part of United States Patent Application Serial 
15 No. 10/728, 136, filed December 4, 2003 which claims the benefit of priority to U.S. 
Provisional Application, Serial No. 60/430,802 filed, December 4, 2002, and is a 
Continuation-In-Part of United States Patent Application Serial No. 10/389,534 filed 
March 14, 2003 which is a Continuation-In-Part of Serial No. 10/389,090 filed March 14, 
2003, which is a Continuation-In-Part of Serial No. 10/143,160 filed May 10, 2002 now 
20 U.S. Patent 6,574,299, which claims the benefit of priority to U.S. Provisional 

Application 60/312,827 filed August 16, 2001, all of which are incorporated by reference 

As previously described, Fig. 1 shows a typical prior art view arrangement of a 
patient on a table that moves through a C-arm device such as the AXIOM Artis MP, 
manufactured by Siemens, that is capable of rotating around the patient, having an x-ray 
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tube source and a detector array, where CB projections are received by the x-ray detector, 
and an image reconstruction process takes place in a computer 4 with a display 6 for 
displaying the reconstructed image. For the subject invention, the detector array is a two- 
dimensional detector array. For example, the array can include two, three or more rows 
5 of plural detectors in each row. If three rows are used with each row having ten detectors, 
then one CB projection set would be thirty individual x-ray detections. 

Alternatively, a conventional gantry, such as ones manufactured by Siemens, 
Toshiba, General Electric, and the like, can be used, as shown by the dotted concentric 
lines, for the X-ray sources and detector array. The gantry can rotate partially, and up to a 
10 full circle, or greater than a full circle. 

Fig. 2 shows an overview of the basic process steps of the invention that occur 
during the image reconstruction process occurring in the computer 4 using a first 
embodiment. An overview of the invention process will now be described. 

A preferred embodiment works with keeping several (approximately 2 to 
1 5 approximately 4) CB projections in computer memory at a time and uses one family of 
lines. 

In the first step 10, a current CB(cone beam) projection set is taken. The next 
steps 20 and 30 identify sets of lines on a virtual x-ray detector array according to the 
novel algorithm, which will be explained later in greater detail. In the given description of 
20 the algorithm the detector array can be considered to be flat, so the selected line can be a 
straight tilted line across the array. 

The next step 40 is the preparation for the filtering step, which includes 
computations of the necessary derivative of the CB projection data for the selected lines. 
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The next step 50 is the convolution of the computed derivative (the processed CB 
data) with a filter along lines from the selected family of lines. This step can also be 
described as shift-invariant filtering of the derivative of the CB projection data. In the 
next step 60, the image of the object being computed is updated by performing back 
5 projection. In the final step 70 it is indicated that the above steps 10-60 are repeated, 

unless image reconstruction is completed or there are no more CB projections to process. 

The invention will now be described in more detail by first describing the main 
inversion formula followed by the novel algorithm of the invention. 

Referring to Fig. 3, we consider a source trajectory consisting of a circle C and a 
10 line segment L attached to C . Fig. 3, shows mathematical notations of the scan. Let 

^ be the point where they intersect. Such a trajectory can be implemented with a C-arm by 

first moving the patient and then rotating the C-arm around the patient. 

We suppose /, 3 5 y(s) e L and 1^3 s-¥ y{s) g C are parameterizations of the 
line and circle, respectively. Here 5* is a real parameter, and y{s) is a point on the 
15 trajectory representing the x-ray source position. We assume that the circle is of radius R 
and centered at the origin. Let U be an open set, such that 
f/c:{(jc,,X2,X3)eC/ +xl <R^}, 

\ 
\ 

Next, we pick a reconstruction point x g t/ , and consider the plane n(x) through 
X and L . n(x) intersects C at two points. One of them is y^ , and the second is denoted 
20 yci^) . Let L^^ix) be the line segment containing x and connecting j^^(x) to L . Then 
yi^{x) e L denotes the other endpoint of L^^{x) . Our procedure determines two 
parametric intervals. The first one /,(jc) c corresponds to the section of L between y^ 
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and y^^ix) . The second one /jCjc) c: /j corresponds to the section of C between and 

We use the following notations in equations 1 and 2 as follows: 
D^iy,®):^ [ f(y-^eOdt,&eS'; (1) 

\x-y(s)\ 

where 

iS^ is the unit sphere, 

/ is the function representing the distribution of the x-ray attenuation coefficient 
1 0 inside the object being scanned, 

0 is a unit vector, 

D^(y,®) is the cone beam transform of / , 

fi{s,x) is the unit vector from the focal point y{s) pointing towards the 
reconstruction point x . 

15 

Next, we suppose first s e Iy{x) . Project x and C onto the detector plane DP{s) 
as shown in Fig. 4. Let h be the height of y{s) above the plane of the circle. We assume 
that DP{s) contains the z -axis and perpendicular to the shortest line segment connecting 
y{s) and DP{s) . It is convenient to introduce the coordinate system on DP{s) , in which 
20 the origin is at (0,0, A) , the first axis is horizontal, and the second axis is vertical. In 
terms of the coordinates on DP{s) that we have just introduced the equation of C is 
given by equation 3. 



' 2 



Consequently, the projected circle C is a parabola, which opens downward. As is seen 
25 from Fig. 4, x. is located above C . By construction, y^^x) e C and is located directly 
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underneath x (see Figures 4 and 5). We find a plane through x and y(s) , which is 
tangent to C at some yX^yX) e I^ix) . This determines the unit vector by equation 4. 

\(y.is,x)-yis))x/3is,x)\' 

5 We now pick a source position y{s) on the circle, j e /2(jf), and define another 

unit vector by equation 5. 

u,(s..):= /.lff-''> x.U.s.I,(x). (5) 

\yis)x/}(s,x)\ 

10 Here y(s) is the velocity vector of the source at the current position. The detector plane 
DP(s) corresponding to the source y(s) e C with various points and lines projected onto 
it is shown in Fig. 6. Using equations (6) and (7) we obtain the following reconstruction 
formula for / gC^((7) by equation 6: 

27r^ k^x I x-y{s) \^ dq ^ ' ^ siny 

15 

where using equation 7 

0* {s, X, := cos rPis, x) + sin ye^ (5, x), 

{s, x) := x) X (5, x). 
X is the cross-product of two vectors, 

q is parameter along the trajectory, 

20 and is defined by equation 8 as follows: 

S^(s,x) = -sgn(M,(5,x) >'(j)), J € /,(x); S^is.x) = l,s e I^(x), (8) 

Suppose, for example, that L is parameterized in such a way that the source moves down 
along I as 5 increases. Then S^{s,x) = 1,5 e /,(x) , If the source moves up along L s 

increases, then (5,x) = -1,5^ e /, (x) . 

25 



12 



C U» S Ql "4-/ A S! Slip 3 & Atty Dkt.: 

UCF-456CIP 

We now describe an efficient (that is, of the convolution-based FBP type) 
implementation of inversion formula equation (6). Pick y(s) e L . For a point x e f/ we 

have to find e /jC^) (see Fig. 5), This determines the filtering line on the detector, 

which is tangent to C at >'(^,) . Here and everywhere below by hat we indicate projection 
of a geometrical object (e.g., circle, line, point, etc.) onto a detector plane. Clearly, all 
other X € C/ which project onto this line to the left of y(Sf) will share it as their filtering 

line. Hence, we can first perform filtering along lines on the detector tangent to C (see 
family in Fig. 7), and then perform backprojection. The range of values, 

s^'" ^Sf< s^°^ , depends on the region of interest (ROI) and is illustrated in Fig. 7. 
Suppose now y{s) e C . 

Under the invention filtering must be performed along lines on the detector 
parallel to y{s) . The resulting family is denoted _2 ^^8- 8. We pick any line from 
_2 • As is seen fi*om Fig. 6, all x whose projection belongs to that line and appears to 

the right of L share it as their filtering line. As before, one can first perform filtering. 
Step 50, (which is shift-invariant) along these lines, and follow it by back projection. Step 
60. 

Let us describe this in more detail. It is clear that s,{s,x) actuedly depends only on 
s and P{s,x) . Therefore, we can write below including equations 9 and 10. 

w* {s, (J, x), {s, py- ( j, /?), k = \,2,/3 = /3(s, x), 

I 



« sm y 



■dy, (9) 
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/W = - AZI , , /'^'f^,, ^.(^>>g(^>^))^^- (10) 

2;r' 15'-'^*<'>|x-j;(5)| 
Any source position y(s) and a filtering line from the corresponding family 
(either or _2 , depending on whether y{s) e L OTy(s) e C ), determine a plane. We call 
it a filtering plane. Since e(Sy /}) - /i = 0j\ e{s^P) |= 1 , we can write equations 1 1 and 12. 
5 (cos e, sin e)\ {s, /?) = (- sin 0, cos <9) (11) 

for all fi,e{s^P) confined to a filtering plane. Here 0 denotes polar angle within a 
filtering plane. Therefore, 



**(^,/?)=f |-^/(j^(^X(cos(^ + r),sin(^ + /))) 



smr 



^ dy (12) 



dq 

for all P confined to a filtering plane. 
10 Equation (12) is of convolution type and one application of Fast Fourier 

Transform (FFT) gives values of ^^{s,P) for all p confined to a filtering plane at once. 

Equations (10) and (12) would represent that the resulting algorithm is of the 
convolution-based FBP type. This means that processing of every CB projection consists 
of two steps. First, shift-invariant and x -independent filtering along a family of lines on 

1 5 the detector is performed. Second, the result is back-projected to update the image 

matrix. A property of the back-projection step is that for any point x on the detector the 
value obtained by filtering at x is used for all points x on the line segment connecting the 
current source position y(s) with x . Since d/dq in (12) is a local operation, each CB 
projection is stored in memory as soon as it has been acquired for a short period of time 

20 for computing this derivative at a few nearby points and is never used later. 
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Now we describe the algorithm in detail following the seven steps 10-70 shown in 
Fig. 2. 

Step 10. We load the current CB(cone beam) projection into computer memory. 
Suppose that the mid point of the CB projections currently stored in memory is yis^) . 
5 The detector plane corresponding to the x-ray source located at yisg) is denoted DP{Sq) . 
Step 20. Fig. 9 is a three substep flow chart for identifying the set of lines, which 
corresponds to step 20 of Fig. 2. Here we assume that the x-ray source is located on the 
line L . Referring to Fig. 7, the set of lines can be selected by the following substeps 2 1 , 
22, and 23. 

10 Step 21 . Choose a discrete set of values of the parameter inside the interval 

s, ,s, ]. 

Step 22. For each 5^, chosen in Step 21 find a line tangent to the projected circle 

c. 

Step 23. The collection of lines constructed in Step 22 is the required set of lines 
1 5 (see Fig. 7 which illustrates the family of lines used in the algorithm of the 

invention). 

Step 30> Here we assume that the x-ray source is located on the circle C . Referring to 
Fig. 8, we form a set of lines parallel to the projection of the circle that cover the 
20 projection of the region of interest (ROI) inside the object being scanned. 

Step 40. Preparation for filtering 
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Fig. 10 is a seven substep flow chart for preparation for filtering, which 
corresponds to step 40 of Fig. 2, which will now be described. 

Step 41. If the x-ray source is located on the line Z , fix a filtering line Ljj^ g_, 

from the set of lines obtained in Step 20. If the x-ray source is located on the 
circle C , fix a filtering line Ljj^ e_2 from the set of lines obtained in Step 30. 

Step 42, Parameterize points on the said line by poleir angle y in the plane 
through yis^) and L^, . 

Step 43. Choose a discrete set of equidistant values /j that will be used later for 
discrete filtering in Step 50. 

Step 44. For each /j find the unit vector /3j which points from yis^) towards the 
point on Z,^, that corresponds to . 

Step 45. Using the CB projection data Dy(y(q),®) for a few values of q close to 
Sq find numerically the derivative (3/ dq)Dj^(y{q),&) \^^^ for all & = /Jj. 
Step 46. Store the computed values of the derivative in computer memory. 
Step 47. Repeat Steps 41-46 for all lines Lj^, . This way we will create the 

processed CB data ^{sq.Pj) corresponding to the x-ray source located at yis^) . 

Step 50, Filtering 

Fig. 1 1 is a seven substep flow chart for filtering, which corresponds to step 50 of 
Fig. 2, which will now be described. 

16 



. p C Ty' U S O '1 e S 3 B Atty Dkt.: 

UCF-456CIP 

Step 5 1 . Fix a filtering line Ljj^ . If the x-ray source is located on the line L , we 

take Ljj, g_ , . If the x-ray source is located on the circle C , we take Ljj^ g_ ^ • 

Step 52. Compute FFT (Fast Fourier Transform) of the values of the said 
processed CB data computed in Step 40 along the said line. 
5 Step 53. Compute FFT of the filter 1 / sin/ 

Step 54. Multiply FFT of the filter 1/sin;^ (the result of Steps 53) and FFT of the 
values of the said processed CB data (the result of Steps 52). 
Step 55. Take the inverse FFT of the result of Step 54. 
Step 56. Store the result of Step 55 in computer memory. 
10 Step 57. Repeat Steps 51-56 for all lines in the said family of lines. This will give 

the filtered CB data ^{s^.Pj) . 

By itself the filtering step can be well known in the field and can be implemented, 
for example, as shown and described in U.S. Patent 5,881,123 to Tam, which is 
incorporated by reference, 

15 

Step 60. Back-projection 

Fig. 12 is an eight substep flow chart for backprojection, which corresponds to 
step 60 of Fig. 2, which will now be described. 

Step 61. Fix a reconstruction point x , which represents a point inside the patient 
20 where it is required to reconstruct the image. 

Step 62. If belongs to I^{x)kjI^{x) , then the said filtered CB data affects the 

image at jc and one performs Steps 63-68. If is not inside /, (x) u {x) , then the 
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said filtered CB data is not used for image reconstruction at jc. In this case go 
back to Step 61 and choose another reconstruction point. 

Step 63. Find the projection x of x onto the detector plane DP(Sq) and the unit 
vector /?(5o,jc), which points from X-^o) towards x. 

Step 64. Identify filtering lines Z^, e_, orL^, e_2 (depending on where the x-ray 

source is located) and points on the said lines that are close to the said 
projections . This will give a few values of <t>(sQ,j8j) for fij close to /}(Sq,x) . 

Step 65. With interpolation estimate the value of <I>(Sq,/3(Sq,x)) from the said 

values of ^{s^^fij) for pj close to /?(5o,jc). 

Step 66. Compute the contribution from the said filtered CB data to the image 
being reconstructed at the point x by multiplying 0{sq^P{Sq,x)) by 

-Sf^{sQ,x)/(27r^ \x-y{sQ) |) . The quantity S^^is^^x) is defined by equation (8). 

Step 67. Add the said contribution to the image being reconstructed at the point x 
according to a pre-selected scheme (for example, the Trapezoidal scheme) for 
approximate evaluation of the integral in equation (10). 
Step 68. Go to Step 61 and choose a different reconstruction point x . 

Step 70, Go to Step 10 (Fig. 2) and load the next CB projection into computer memory. 
The image can be displayed at all reconstruction points x for which the image 
reconstruction process has been completed (that is, all the subsequent CB projections are 
not needed for reconstructing the image at those points). Discard from the computer 
memory all the CB projections that are not needed for image reconstruction at points 
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where the image reconstruction process has not completed. The algorithm concludes 
when the scan is finished or the image reconstmction process has completed at all the 
required points. 

The invention is not limited to an object that undergoes a scan consisting of a 
5 single circle and a single line segment. The algorithm can be applied to trajectories 
consisting of several circles and line segments by applying it to various circle and line 
pairs, and then combining the results. The algorithm can be applied to trajectories in 
which a planar curve is not necessarily a circle, but, for example, an ellipse, and the like. 

10 Other Embodiments of the invention are possible. For example, one can integrate 

by parts in equation (6) as described in the inventor's previous U.S. Patent Application 
Serial No. 10/143,160 filed May 10, 2002 now U.S. Patent 6,574,299, now incorporated 
by reference, to get an exact convolution-based FBP-type inversion formula which 
requires keeping only one CB projection in computer memory. The algorithmic 

15 implementation of this alternative embodiment can be similar to and include the 

algorithmic implementation of Embodiment Two in the inventor's previous U.S. Patent 
Application Serial No. 10/143,160 filed May 10, 2002 now U.S. Patent 6,574,299, now 
incorporated by reference. The corresponding equations v^ll now be described. 

20 Introduce the following notations in equation 1 3 : 
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While the invention has been described with rotating C-arm type devices, the 
invention can be used with rotating gantry devices. 

Furthermore, the amount of rotating can include a single rotational curve of at 
least approximately 5 degrees up to approximately 360 degrees or greater. Theoretically, 
there is no limit on the minimum range of rotation. Under realistic practical 
circumstances, a minimum range of rotation is between approximately 10 and 
approximately 20 degrees. 

The circle and line scanning of an object can have a line scanning before or after a 
single rotational curve scan as defined above. 

Subsequent circle and line scanning can occur as needed for image reconstruction. 

Although the preferred embodiments describe applications of using x-ray sources 
for creating data for image reconstruction, the invention can be applicable with other 
sources that create line integral data for image reconstruction, such as but not limited to 
early arriving photons. 

While the invention has been described, disclosed, illustrated and shown in 
various terms of certain embodiments or modifications which it has presumed in practice, 
the scope of the invention is not intended to be, nor should it be deemed to be, limited 
thereby and such other modifications or embodiments as may be suggested by the 
teachings herein are particularly reserved especially as they fall within the breadth and 
scope of the claims here appended. 
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